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Using a relativistic Dirac-Brueckner analysis the pion contribution to the ground state energy of 
| nuclear matter is studied. Evidence is presented that the role of the tensor force in the saturation 

. mechanism is substantially reduced compared to its dominant role in a usual non-relativistic treat- 

Q\ ' ment. The reduction of the pion contribution in nuclear matter is due to many-body effects present 

in a relativistic treatment. In particular, we show that the damping of OPEP is actually due to the 
decrease of M* /M with increasing density. 
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Within a non-relativistic framework the tensor force plays a dominant role both in the formation of the deuteron 
and in nuclear matter saturation mechanism. For the deuteron 



( i D 1 | T/ total | 3 Si) ~ 2{ 6 D X | y tensor | 3 Si) ~ -24MeV, 

while the kinetic energy is ~ 22 MeV, adding up to a binding energy of ~ — 2McV. In a non-relativistic Bethe-Brucckner 
^ , calculation of nuclear matter one finds typically [|l| 

(NM | K I NM) non _ rclativistic ~ -34 (p/ Po ) OA5 MeV, (1) 

0^ 1 where p = 0.17 fm 3 is the saturation density of normal nuclear matter. The exponent of p is markedly less 
than the nominally expected value of 1 because of Pauli blocking. It is more than four decades ago that it was 
pointed out that in the absence of the tensor force the saturation properties of normal nuclear matter cannot be 
understood assuming the validity of a non-relativistic quantum description. In particular, saturation occurs at a 
density one order of magnitude or more higher than po- The effect of the pion, being given essentially by the second 
order in the one-pion-exchange potential (OPEP) contribution, is attractive and large at low density. Pauli blocking 
reduces the attraction as density increases. This mechanism leads to saturation of nuclear matter at a density well in 
. j-h , correspondence with the empirical data. The above conclusions are obviously based on a non-relativistic description 
1 of the nuclear many body system. 

The relativistic results for the contribution of V„ to the deuteron || is (D | V w | D) = — 31MeV, suggesting an equally 
important role of the pion. In sharp contrast, this seems not to be the case in a relativistic treatment of nuclear matter. 
Strong scalar (S) and vector (V) fields^] of the order of a few hundred MeV are typical for relativistic theories || ||] 
based on a meson theoretical description of the nuclear force. These values are consistent with expectations based 
on the studies of scattering of ~ 1 GeV protons by nuclei. The large scalar fields have far reaching consequences in 
nuclear matter through the strongly medium modified nucleon mass M* = M + S. The saturation mechanism is 
believed to rest upon the decrease of M* with increasing density. It is the only possible mechanism for saturation 
m a mean field theory (MFT) like the QHD @. A Dirac-Brueckner (D-B) analysis J|| is at present the best 
tool we have for a relativistic study of nuclear matter. In Ref. || one finds Sd-b = —306 (p/po) ' 81 MeV, while 
Vd-b — 233 (p/po) ' 97 MeV. This reduction in the rate of increase of the strength of S compared to the increase of V 
with increasing density is the saturation mechanism in D-B. In a MFT the vector field and the scalar field arise from 
the exchange of the u> and the a meson. A MFT treatment using the interaction of Ref. ||, with only the a— and 
w-exchange contributing, yields Smft ~ — 358 {p/po) ' 92 MeV and Vmft ~ 295 MeV. The reduced rate of increase 
of Smft with increasing p is entirely due to the decrease of scalar charge of the nucleon. In a D-B study, where 



lr To be precise, here we refer to the scalar and vector parts of the self-energy as fields. 
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ladders of meson exchanges are summed, there are other factors contributing to both S and V. As a result there 
are significant differences in both magnitude and density dependence of S and V in the two treatments. Thus when 
ladders are summed the saturation mechanism need not be exclusively due to er-exchange. 

Because of the dominant role of OPEP in the saturation mechanism in non-relativistic studies of nuclear matter, 
one should expect a similar role in relativistic studies. The important role of the density dependence of S does not 
exclude it. In this letter we examine the role of OPEP in D-B and show that it is substantially reduced due to 
relativity. Since the contribution of OPEP to the deuteron binding energy remains large in a relativistic treatment 
the damping in nuclear matter must be due to many-body effects. We find that it can be attributed to the decrease 
of M* /M with increasing density. Assuming that the nuclear matter is uniform in space and constant in time with 
a given density p, the baryon current is given by — pu^ with it M = (0, 1) being the unit vector in the nuclear 
matter frame. Relativistic covariance implies that the self-energy contribution of the nucleon with momentum p can 
be characterized by 

E = S s - S" 7 • u - E" 7 • p. (2) 
The medium modified mass of the nucleon is then given by 

M* = (M + E S )/(1 + E W ). 

Scaling all relevant momenta in terms of M* units we may express the ground state binding energy E/A in terms of 
a dimensionless NN G-matrix G as 

o 2 f k F /M* ,3f 

+ / r%^3 / 7^^— r-(^i,Ai,ii;^2,A 2 ,i2 | G | £i,Ai,ii;£ 2 ,A 2 ,i 2 )], (3) 

where £ n — p n /M* are the scaled momenta of the two interacting nucleons and E(l) = \Jl 2 + 1. The G-matrix satisfies 
the Dirac-Bethe-Brueckner-Goldstone equation involving dimensionless quantities. Within a relativistic quasipotcntial 
approach it has in the NN cm. system the form 

(i 1 1 e) = {? I v I i) 
dH" 



+ E I $\V\V) f\G\l), (4) 

where we sum over the p-spin, helicities and isospin and integrate over the scaled relative momenta I" of allowed 
intermediate states. Furthermore, S% is the 2-nucleon Green function, including the Pauli- blocking operator Qp au u- 
Using a Blankenbecler-Sugar-Logunov-Thavkhclidzc prescription §2 has the form 



+ 

[E no - fc-f + M*] [Efjo + fc-7 + M *] 
{E*+E f f{El-E f -ie) 



Qpauli (5) 



with k = l"M*, Ef = Wq/(1 + E 11 ) and Wq — (pi +P2) 2 , Wo being the total invariant energy of the final state. 
Furthermore, in Eq. (Q) V are the dimensionless quasi-potential matrix elements in Dirac space 

(? I V\ £) = M* 2 

where it's are the positive and negative energy (p = ±1/2) spinors for mass M* fermions, satisfying 

( P E* lQ - 7 ■ p - M*)u XjP (p) = 0. (7) 

with E* = (p 2 + M* 2 ) 1 / 2 . Interdependence of the nucleon self-energy and the G-matrix, for a fixed given nuclear 
matter density, requires selfconsistency in solving the D-B equations. The NN quasipotcntial is taken to be given 
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by the relativistic one-boson-exchange (OBE) model with tt, p, us, a, 6 and ?7-mesons. We have, ignoring the isospin 
factors, for the scalar, vector and pseudo-scalar meson exchanges in nuclear matter 

(^ I V u , p I p) = ^[n(^){ 7/J + i^^ip' ~ PY}U(P)] [1) ^ ^ + ^ ~ ^<^V ~ P)M~P)\ W , (9) 

{0 \v P \p} = ~(^y{u(jf) l5 (f-ftu(p)}W ^_^ 2 + m2p H~i?h 5 (f-p'M-p)}W. (io) 

We will assume, that possible density dependence of meson-nucleon couplings and meson masses can be neglected. 
Using an angular-averaged Pauli-blocking operator we may write for the 2-nucleon Green function 

M* 2 1 , . 

(11) 



E* 2 (P,p) - E* 2 (P,k) (p/M*) 2 - (k/M*) 2 

The last line of Eq. ( pi] ) indicates that the range of the intermediate momentum k is controlled by M* . The smaller 
the M* , the smaller is the contributing range of k. This fact has an important effect on the n and r/ and the Pauli 
coupling of the vector mesons. The corresponding matrix elements have two extra powers of k in the numerator 
compared to other matrix elements. So with increasing density and, therefore, with decreasing M* these matrix 
elements have suppression factors not present in all other matrix elements. 

Given the solutions of the G-matrix the binding energy is calculated using Eq. (|3|). The potential contributions of 
the pion can be calculated using the Hellmann-Feynman theorem 

(NM | V4 | NM) = gl NN —l—(E/A), (12) 

9ttN N 

and we find, that the contribution of to E/A is 

(NM | | NM) Rc i ativistic ~ -20 (p/po)°- 16 MeV. (13) 
It is considerably suppressed compared to the value given by Eq. (0) for the non-relativistic case. 
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FIG. 1. Plots of the quantities E/A (curve a) and E/A - (NM | w w | NM) ~ 17 MeV (curve b). The numbers are based on 
the work of Ref. |j. 

Furthermore, it has only a minor role in the saturation mechanism. This is exhibited in Fig. [l] where we plot E/A 
(curve a) and E/A — (NM \ v v \ NM) — 17 MeV (curve b), all obtained from ||. The two curves have practically 
the same density dependence verifying that OPEP contributes little to the saturation mechanism. The subtraction 
of 17Mev in curve (b) makes the scale more compact. 
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The above results can be understood qualitatively by examining the 2nd order contributions to the G-matrix. 
Keeping only the positive energy M* state contributions in the intermediate states we have 



*(t'\V I V) . ™ T (V' | V 1 1), (14) 
Wo - L 2 /4 - I" 2 

where Wq = Wq/M* and L = (pi +f>2)/M* . Now the tensor force contributes only to the second term of Eq. (pl|). 
It has the form 

(? I % I I) = -(^) 2 (^)WW ~f)u(M^ 

*ar^L — o fcM~W -f)<-?)} {2 \ (15) 

(£' - £) 2 + ml 

where u's are the positive energy spinors. Similarly as in the non-relativistic case the 2nd order pion contribution is 
strongly density dependent because of the Pauli blocking. However as is exhibited in Eq. (jl^) in the relativistic case 
the effective coupling of the pion to the nucleon is suppressed by a factor of 

The ^jt- suppresssion can be corroborated in more detail by the following calculation. Let us modify the S from 
the self-consistent calculation of Ref. || by multiplying it with the factor a < 1 thus generating a M* = M + aS. 
By using the modified scalar self-energy in the nucleon propagators we recalculate first the G matrices and then E/A 
and finally (NM | V n \ NM) using Eq. (|l2|). Only the a — 1 analysis is self-consistent, others are not. But such 
a calculation is particularly suitable to exhibit the role of M* /M on the OPEP contribution. The Fig. ^ exhibits 
clearly the damping due to decreasing M* /M. We stress that the mechanism of damping is generic to any relativistic 
treatment and not particular to either Ref. || or the use of Ref. |l(J . 



-15 




-45 1 

0.08 0.13 0.18 0.23 

FIG. 2. Plots of {NM | Vn \ NM) with the parameters of Ref. §). In the G-matrix calculations S is replaced with aS. The 
plots are for a = 0., 0.5 and 1.0. The last one is the result of self-consistent calculations of Ref. Q|. The other two are not 
self-consistent. 



We want to be careful that the present work not be interpreted as providing support for MFT. As shown in Fig. |3|, 
the actual strengths of the scalar and vector potentials found in the OBE interaction does lead in a MFT treatment to 
distinctly different predictions from D-B calculations for E/A. Undoubtedly, if one releases oneself from the constraint 
of fitting NN data and freely chooses the NN interaction one can obtain proper binding and saturation of nuclear 
matter with a MFT calculation. 
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FIG. 3. Plots of E/A from the D-B calculations of Ref. || and from a MFT calculation with the same quasipotential. 

In conclusion, we have established that in the relativistic treatment the tensor force contributions are reduced in 
size in nuclear matter. Because of this, in complete contrast to the non-relativistic situation, they cease to play a 
role in the saturation mechanism. The reduction of the tensor force contributions is principally due to the relativistic 
M*/M effect, as illustrated with Fig. ||. But even the reduced role of OPEP is not negligible. As noted, it contributes 
— 20MeV to E/A. This and an enormous body of pre-existing evidence lead us to conclude that there is no viable 
non-relativistic theory of nuclear matter based on the use of NN forces which fit two-body data which is compatible 
with the relativistic theory. The dominant mechanism of saturation of nuclear matter is basically very different in the 
two approaches. In the non-relativistic approach it is the density-dependent reduction due to Pauli blocking of the 
attraction from tensor force, while in the relativistic approach it is the reduction of the rate of growth with increasing 
p of the attraction from the scalar field relative to the growth of repulsion from the vector field. 
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